Enumeration of 4 x 4 magic squares

A magic square is an n × n array of distinct positive integers whose sum along any row, column, or main diagonal is the same number. We compute the number of such squares for n = 4, as a function of either the magic sum or an upper bound on the entries. The previous record for both functions was the n = 3 case. Our methods are based on inside-out polytopes, i.e., the combination of hyperplane arrangements and Ehrhart’s theory of lattice-point enumeration. 1. It’s a kind of magic A magic square is an n × n array of distinct positive integers whose sum along any row, column, or main diagonal is the same number, the magic sum. The history of magic squares is well documented; see, e.g., [8, 9, 21]. The contents of a magic square have varied with time and writer; usually they have been the first n consecutive positive integers, but often any arithmetic sequence and sometimes fairly arbitrary numbers. The fixed ideas are that they are integers, positive, and distinct. In the last century mathematicians took an interest in results about the number of squares with a fixed magic sum, but with simplifications: diagonal sums were often omitted and the fundamental requirement of distinctness was almost invariably neglected [1, 3, 11, 16]. For example, classical formulas of MacMahon [15] include